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Abstract 

This paper is concerned with the topological invariant of a graph given by the maximum 
degree of a Markov basis element for the corresponding graph model for binary contingency 
tables. We describe a degree four Markov basis for the model when the underlying graph 
is a cycle and generalize this result to the complete bipartite graph K2.n- We also give a 
combinatorial classification of degree two and three Markov basis moves as well as a Buchberger- 
free algorithm to compute moves of arbitrary given degree. Finally, we compute the algebraic 
degree of the model when the underlying graph is a forest. 

Keywords: Markov bases, contingency tables, graphical models, hierarchical models, toric 
ideals. 

1 Introduction and Definitions 

The study of multidimensional tables and their marginals is of central importance whenever one 
wishes to make inferences based on statistical samples. In general, one is presented with a nonneg- 
ative integral table of data of size di x ■ ■ ■ x dn and a simplicial complex A on { 1 , . . . , n} which 
encodes the specific marginals we would like to compute; this is called a hierarchical model. 

Certainly, the oldest example of such a model is the case of computing the row and column 
sums of a matrix. If the matrix is a square m x m matrix and we require that all the row and 
column sums have the same value k, one calls such a matrix a semi-magic square with magic sum 
k. Here, our simplicial complex consists of two isolated points. 

In the more general statistical situations, each node of the simplicial complex corresponds to a 
feature of a population sample (e.g. eye color) and the levels of the table correspond to different 
states of the feature (e.g. green, brown, hazel, etc.). The faces of the simplicial complex are intended 
to model interactions between the features. One of the most fundamental questions in statistical 
analysis is: do the data appear to be satisfied by a given model? One way to test the hypothesis is 
to compare the sample data to the maximum likelihood estimate using the or statistic. One 
problem with applying this approach directly is that the data is always an integral table, while the 
maximum likelihood estimate almost never is. As a result, there might be no integral table with the 
same marginals which has a small statistic. This problem is especially dramatic when analyzing 
the large, sparse data sets which occur in real world situations (e.g. census data). To remedy this 
situation, one may attempt to decide whether or not a model fits the data by comparing statistics 
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of the data table with statistics of random integral tables with the same marginals. If the statistic 
of the table of data is exceptional one could hope to conclude that the table was also exceptional. 
For example, if the statistic of the table of data was exceptionally small, one could conclude 
that the data did, in fact, fit the model. 

We are now left with the problem of generating random integral tables from the set of all 
nonnegative integral tables with fixed marginals. One solution is to perform a random walk over 
the set of all nonnegative integral tables with given fixed marginal. Such a random walk can be 
taken by first finding a suitable set of "moves" (these are tables with integral entries which have all 
their marginals equal to zero) and randomly adding moves to some starting table. It is at this point 
in the story that computational commutative algebra enters the picture: finding such a set of moves 
is equivalent to finding a generating set for the associated toric ideal. For a detailed introduction 
to the connections between toric algebra and multidimensional contingency tables, see [2], 0) and 

EI- 

In this paper, we restrict attention to tables where di = 2 for all i and for which the underlying 
simplicial complex is a graph; that is, we compute only two- and one-way marginals of our binary 
table. We refer to such models as binary graph models. These are generally not the usual graphical 
models studied so frequently in statistics, where the simplicial complex consists of the cliques of 
the underlying graph G ^01 • Our notion of graph model coincides with the more familiar graphical 
model if and only if the graph has no three-cycle. Now, we will give two formal presentations of 
the objects of interest in this paper. 

Let G be a graph on the n-element vertex set [n] := {1, 2, . . . , n} with edge set E{G). Denote 
by Iso{G) the set of isolated vertices of G. For each edge {j, k} of the graph G consider the linear 
transformations TTj^fc 

TTj^fc : — > 

and for each isolated vertex k of G consider the linear transformations nk 

TTk : Z^" — > I? 

Gil. ..in ' ^ ^ik ■ 

We think of the maps iTj^k as computing the 2-way marginal of a 2 x • • • x 2 table corresponding to 
the edge {j, k} and the maps vr^ as computing the 1-way marginal of a 2 x • • • x 2 table corresponding 
to the vertex k. We define the map ttq by taking all the marginal computations induced by a given 
graph as 

{j,k}eE{G) k£lso{G) 

V' > 7rj-fc(t>) TTkiv). 

{j,k}<^E{G) kelso(G) 

We say that ttg is the map which computes the marginals of a 2 x • • • x 2 table according 
to the graph G. The matrix which represents this linear transformation will be denoted Aq and 
the polytope which is the convex hull of the columns of Aq is denoted Pq where we consider the 
columns of Aq as vectors in M*^ for an appropriate d. A move for G is an element of the integral 
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kernel of ttg] that is, a move is an integral table which does not change the G-marginals of a table 
it is added to. In general we are interested in sets of moves with special properties. 

Definition 1.1. A finite subset of moves B C kerx(7rG') is called a Markov basis for the graph 
G if for every pair of nonnegative integral tables vi,V2 G N^" with the same G-marginals ircivi) = 
ttg{v2), there is a sequence of moves {ui}[^i C =bi? such that 

I 

Vl + y^^Ui = V2 

1=1 

and 

j 

vi + J2'^i^ N^" for all l<j<l. 

i=l 

There is also a much shorter and more algebraic way to arrive at this definition. Recall that G 
is a graph on the n-element vertex set [n] := {1, 2, . . . , n} with edge set E(G), and isolated vertices 
Iso{G). Consider the map of polynomial rings 

: CK,....Ji, G {0, 1}] C[tg£,tSf |{i, A;} G E{G) and / G Iso{G)] 

Ph.. .In II J_| • 

{j,fc}G£;{G) leIso{G) 

The object of interest in this paper is the ideal denoted Iq = ^er{<j)G), which we call the ideal of a 
binary graph model. It is a toric ideal: a prime ideal generated by monomial differences whose 
leading and trailing terms have disjoint support. 

Markov bases and toric ideals are connected by the following fundamental theorem. 

Theorem 1.2. A finite subset of moves B = {nj}'^| C kerz(7rG) is a Markov basis for G if 
and only if the set of binomials {p^i — p'^i } is a generating set for Iq- 

Here, we write Ui = uf — u\ as the difference of two positive vectors of disjoint support. In 
light of Theorem 11.21 we will use the expressions "Markov basis for G" and "generating set for 
/g" interchangeably throughout this paper. Similarly, we can interchange the words "move" and 
"binomial" whenever we are discussing the Markov bases/generating sets of Ig- These definitions 
are best illustrated by a simple example. 

Example 1.3. Consider the graph G on four nodes with two edges {1,2} and {2,3} and one 
isolated vertex 4. The map (j)G is a map from a polynomial ring in sixteen variables to a polynomial 
ring in ten variables. It is the map of rings 

(/•G : 'C[pijki\i,j,k,l G {0,1}] — > C[rij,Sjk,ti] 

Pijkl ' *■ ^ij ' Sjk ' tl- 

On the other hand, the marginal map ttg is given by a 10 x 16 matrix Ag- It is the matrix 
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A 



G 



/llllOOOOOOOOOOOON^ 
0000111100000000 
0000000011110000 
0000000000001111 
1100000011000000 
0011000000110000 
0000110000001100 
0000001100000011 
1010101010101010 

Voioioioioioioioi/ 

The polytope Pq in M}^ has dimension 6 and has 10 facets. These facets are indexed naturally 
by the rows of Aq and the facet defining inequahties are given by > with one inequality for 
each row. The ideal Iq has a Markov basis consisting of quadratic moves. These are 



PojOhPijih - PojihPijOh with j, /i,/2 G {0,1} 



and 



PhjikihPi2j2k2l2 ~ Phjlkil2Pi2j2k2li 'with il,jl,ki,li,i2,j2,k2,l2 G {0,1}. 

These generators are also a Grobner basis with respect to the reverse lexicographic term order with 
po-o < < Pl-l- 

In general, when the underlying graph is a forest, the toric ideal Ig is relatively well understood. 

Theorem 1.4. [H MSf The ideal Ig is minimally generated by quadrics if and only if G is a 
forest. In this case, the set of quadratic squarefree binomials in Ig forms a Grobner basis with 
respect to the reverse lexicographic term order with po-- o Po-- oi -<■■■-< Pi --i- 

In this paper, we are primarily concerned with investigating graphs which contain cycles. One 
fundamental question is to compute the following invariant of a graph. 

Definition 1.5. Let G be a graph. The Markov width is the degree of the largest minimal 

generator of the toric ideal Ig- 

Studying the Markov width of a graph is of fundamental importance for statistical applications 
because it relates the complexity of analyzing data to the complexity of the underlying graphical 
structure. Note that. Theorem ^31 states that the graphs with fJ.{G) = 2 are precisely forests, which 
are certainly topologically simple. Indeed, the Markov width of a graph G is topological in nature, 
by which we mean that /u(G) can only decrease under the operations of vertex deletion and edge 
contraction; we will show this in Section 4. Since these operations interact nicely with the toric 
ideals of the initial and final graph, we will use the following definition throughout. 

Definition 1.6. Let G be a graph. By a minor ofG, we mean a graph H which can be obtained 
from G via a sequence of edge contractions and vertex deletions. 

This is different from the usual definition of a graph minor in that we do not allow edge deletion, 
whose interaction with the toric ideal is more complicated. 
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The rest of the paper is organized as follows. In the next section, we discuss computational 
results for graphs with few vertices and suggest some conjectures based upon these data. In 
particular, we have computed the Markov width for all graphs on five vertices and many of 

the graphs on six vertices. In the third section, we prove that the n-cycle and K2^n have Markov 
width 4. Furthermore, we are able to explicitly describe the moves needed in the Markov bases for 
these graph. The fourth section is devoted to the inverse problem: that is, studying which graphs 
may have Markov basis elements of a given degree. We give an algorithm which does not depend 
on computing S'-pairs for computing all the minimal generators of a given degree for the ideals Iq- 
As a consequence of our algorithm, we give a combinatorial characterization of moves of degree two 
and three. In the final section, we return to the study of forests and use the reverse lexicographic 
Grobner basis from Theorem 11.41 to derive combinatorial formulae for the algebraic degree of Iq 
whenever G is a forest. 

2 Graphs with Few Vertices 

In this section we discuss and display computational results about the ideals of binary graph 
models. In particular, we describe generating sets for the ideals Iq for all graphs G with fewer 
than five vertices and all the graphs on six vertices with at most eight edges. These computational 
results suggest many conjectures and open problems which we describe at the end of the section. 
All of our computations were carried out using the toric Grobner basis program 4ti2 [Hj and the 
computational algebra system Macaulay 2 [Jj. We limit our description to graphs which cannot 
be "glued" together from smaller graphs based on the following definition and theorem. 

Definition 2.1. Let (Vi, 5, V2) he a partition of the vertex set of a graph G such that 

1. there are no edges in G between Vi and V2 <ind 

2. S is either the empty set or S is a common vertex or edge of the induced subgraphs Gi and 
G2 with vertex sets U 5 and V2 U S, respectively. 

Then G is called reducible with decomposition {Vi, S, V2). 

Theorem 2.2. Q \^ Let G be a graph which is reducible and Gi and G2 the induced subgraphs 
arising from the vertex decomposition. Then there is a degree preserving operation with which one 
can build generating sets and Grobner bases for Iq from generating sets and Grobner bases of Id 
and Ig2- In particular, fJ-{G) = max(/i(Gi), ^(^2)). 

There are precisely one graph on three vertices, two graphs on four vertices, six graphs on five 
vertices, and six graphs with six vertices and at most eight edges which are not reducible. We will 
briefly describe these graphs and their Markov bases. 

2.1 Three and Four Vertices 

The only graph on three vertices which is not reducible is the complete graph K^. The Markov 
basis of IK3 consists of the single degree four binomial 

PoooPoiiPioiPiio - PooiPoioPiooPiii- 
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The two irreducible graphs on four vertices are C4 and i^4. The graph C4 is the four cycle with 
E{Ci) = {12, 23, 34, 14}. The Markov basis of C4 consists of eight quadrics such as 

PooooPoioi — PoooiPoioo 

and eight quartics. The complete graph i^4 has Markov basis consisting of 20 moves of degree four 
and 40 sextic binomials such as 

2 2 
Poooo^'oiiiPioiiyiioiPiiio — PoooiPooioPoiooPioooPiiii- 

2.2 Five Vertices 

There are six graphs on five vertices which cannot be decomposed into subgraphs. These are 
the graphs we denote C5, i^2,35 SP, BP, and K^. The graph C5 is the five-cycle E{C^) = 
{12,23,34,45,15}. Its Markov basis consists of 80 quadrics and 40 quartics. The graph K2,z 
is the complete bipartite graph E{K2fi) = {13,14,15,23,24,25}. Its Markov basis consists of 
44 quadrics and 420 quartics. The graph i^'4 is the graph obtained from K/^^ by subdividing an 
edge, E{Ki) = {12,15,23,24,34,35,45}. The Markov basis for I^^ consists of 32 quadrics, 473 
quartics, and 160 sextics. The graph SP is the edge graph of the square pyramid, E{SP) = 
{12, 13, 15, 23, 24, 34, 35, 45}. The Markov basis of SP consists of 16 quadrics, 671 quartics, and 320 
sextics. The graph BP is the edge graph of the bipyramid over a triangle, E[BP) = {12, 13, 15, 23, 
24,25,34,35,45}. Its Markov basis consists of 8 quadrics, 436 quartics, and 2872 sextic binomials. 
Finally, is the complete graph on five vertices. The Markov basis of consists of 260 degree 
four moves, 3952 sextic binomials, 846 binomials of degree eight such as 

3 3 

PoOOOO^'OllllPlOlllPllOllPlllOl^'llllO — POOOOlPOOOlOPOOlOO^^OlOOOPlOOOOPlllll 

and 480 degree ten binomials like 

2222 4 
Pooooo^'oiiii^'ioooiJ'iooioPioiooPiiooo - PoooioPooioiPoiooiPoiiioPioooo^'iooiiPiiiii • 

2.3 Six Vertices 

There are a total of 29 graphs on six vertices which are not reducible. We were able to compute 
Markov bases for the six irreducible graphs on six vertices which have at most eight edges. It 
remains a major computational challenge to determine Markov bases of the other 23 irreducible 
graphs on six vertices. The six irreducible graphs on six vertices with less than nine edges will be 
denoted Cg, -fC2,4, G151, 6*153, and G154. 

The graph Cg is the six cycle with edge set E{Cq) = {12, 23, 34, 45, 56, 16}. The graph i^2,4 is 
the complete bipartite graph with edge set E{K2^a) = {13, 14, 15, 16, 23, 24, 25, 26}. The remaining 
graphs do not have special names: the labels we have chosen come from |T2- These four graphs have 
edge sets E{Gi29) = {12, 15, 23, 26, 34, 45, 56}, E{Gi5i) = {12, 14, 23, 26, 34, 36, 45, 46}, ^(Gi53) = 
{12, 15, 16, 23, 24, 45, 46, 56}, and E{Gir>i) = {12, 14, 23, 25, 34, 36, 45, 56}. The data regarding the 
Markov bases of these graphs as well as all the irreducible graphs on five and fewer vertices is 
summarized in the following table. The columns are labeled by the particular irreducible graph, 
the rows are labelled by degree of minimal generators and the table entries are the number of 
minimal generators of a given degree. 
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K3 


Ca 


K4 


C5 


K2,3 




SP 


BP 


K5 


Ce 


K2,A 


G129 


Gi5i 


^153 


^154 


2 





8 





80 


44 


32 


16 


8 





528 


236 


360 


280 


320 


256 


4 


1 


8 


20 


40 


420 


473 


671 


436 


260 


160 


11696 


2636 


4949 


4149 


7784 


6 








40 








160 


320 


2872 


3952 











640 


480 


640 


8 


























846 




















10 


























480 




















Total 


1 


16 


60 


120 


464 


665 


1007 


3316 


5538 


688 


11932 


2996 


5869 


4949 


8680 


^(G) 


4 


4 


6 


4 


4 


6 


6 


6 


10 


4 


4 


4 


6 


6 


6 



For these graphs, ah generators are in even degree. This is not true in general, however, as we 
will demonstrate in Section 4. Theorem II .41 characterizes graphs with = 2 as forests, but the 
next case is already quite interesting. 



Problem 2.3. Characterize those graphs with Markov width fJ,{G) = 4. 

In the next section we will show that cycles and the complete bipartite graphs K2^n have Markov 
bases consisting of moves of degree four or less, but from the data we see that this is not yet a 
complete characterization. 

A natural class of graphs which one would hope to understand is planar graphs. The data above 
suggest the following optimistic conjecture. 

Conjecture 2.4. There is a universal constant C such that the Markov width fj,{G) < C whenever 
G is a planar graph. Even stronger, C = 6. 

On the other hand, the data also suggest the following conjecture. 

Conjecture 2.5. The invariant fJ,{G) is a function only of the tree width of G. 

The tree width is a topological invariant of a graph G which is equal to one less than the size 
of the largest clique in the chordal graph containing G which has the smallest maximal clique. For 
example, forests are precisely those graphs with tree width zero or one, and indeed Theorem 11.41 
tells us that these graphs all have Markov width two. Conjecture 12. 51 also agrees with Theorem 12. 2| 
since the tree width of a reducible graph is the maximum of the tree widths of its components. 

While the limited information we have suggests both Conjecture 12.41 and Conjecture 12. 5| they 
cannot both be true: there are planar graphs with arbitrarily large tree width. For example, grid 
graphs can have arbitrarily large clique size in their minimal chordal triangulations. This suggests 
another research problem. 

Problem 2.6. Study the binary graph model Iq for the family of m x n grid graphs. 

We do know that can be arbitrarily large. For example, for the complete graph Km we 

can construct generators of large degree. 

Proposition 2.7. The complete graph Km with m > 3 has Markov width fj,{Km) > 2m — 2. 

Proof. It suffices to show that there is a minimal generator of Km of degree 2m — 2. For this 
consider the binomial 
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i=l i=l 

where is the string of all zeros, 1 is the string of all ones, and Cj is the ith unit vector. Then 
the monomials coming from the leading and trailing terms are the only monomials which have the 
given image under . Equivalently, the corresponding tables are the only two tables which have 
these same fixed marginals under vrx^- Since the leading and trailing terms have disjoint support, 
this binomial must appear in every Markov basis for Ik^- ^ 

Of course, this bound is already not tight for m = 5, where it yields niK^) > 8 despite the fact 
that has Markov width 10. In general we suspect that ^{Km) grows exponentially in m. 

3 Cycles and Bipartite Graphs 

In this section we confirm the observations from the second section: the ideal of the cycle and the 
complete bipartite graph -fC2,n s^re generated in degrees two and four. 

3.1 Cycles 

For ease of notation, we will represent a binomial such as pioiiPiiio— PiiiiPioio in tableau notation 
as 



10 11 
1110 



1111 
10 10 



The tableau are obtained from a binomial by recording the indices of each variable which appears 
in the monomial, repeating indices when a variable appears to a power greater than one. We say 
that one binomial contains another if it does so in the Graver sense; that is, p" — contains 

— p** if a < u and b < v componentwise. 

We first prove the following theorem bounding the degree of minimal generators for the n-cycle. 

Theorem 3.1. Let Cn he the n-cycle graph. Then fJ-{Cn) = 4, and in particular Ic„ is generated 

in degrees 2 and 4- 

Proof. We will start with an arbitrary binomial / in the ideal, and express it as a linear combination 
of elements either of lower degree or of degree at most 4. 

Given any binomial, take one variable from each monomial such that the two variables chosen 
agree in first and last index. Our strategy will be as follows: by adding multiples of ideal elements 
of degree 4 and less, we will eventually obtain a binomial in which both monomials have the same 
variable. Dividing out by this variable (which clearly does not affect membership in Iq) yields a 
binomial of lower degree which must still be in the ideal, completing the proof. 

We now start this process. We have a variable pi?...?i in the first term of /, and a variable 
in the second term. We wish to eliminate all disagreements between these indices by "moving" the 
table entries corresponding to these binary strings using binomials of degree 4 or less. 

Consider any block of disagreements, in which, without loss of generality, the indices of these 
variables look like (• • -TIO • ■ • 01? • • • ) and (• • -Til • • • 11? • • • ). We propose to add some multiple of 
an ideal element of degree g at most 4 so that the resulting binomial contains the two variables in 
question, and is unchanged except that some of the disagreements in the block have been removed. 
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The two sets of index strings in g will agree on the portion of C„ outside of the block in question, 
counting the boundary elements. Essentially, we are performing a local move by changing indices 
on a subgraph of C„. Let / represent the image of / in this subgraph, i.e. under the map sending 
p. ..717... to Pi, where I is the index substring on the block we have. 

Continuing in this manner, by induction on the number of disagreements we eventually obtain 
a binomial for which the same variable appears in both monomials, completing the proof. We now 
construct the ideal element g which we will add a multiple of. We first construct the part g which 
corresponds to the block in question; in the tableaux that follow, we consider only the indices 
corresponding to this block. 

Because this element is in the ideal Ic„, considering the marginal in the first two directions, 
since an element in the first term of / has a 10 marginal, so must an element in the second term. 
So / contains 



[ 1 ••• 1] 



1 1 1 
1 ? 



1 1 

? ? 



Now, if any of the unspecified elements is 1, we let g be the binomial which switches the 
intervening substrings, i.e. something of the form 



11111 
1 1 



10 11 
1 1 1 1 



where A is the remainder of the index string of the element that the second term of / contains 
starting with 10. We fill in the rest of both terms of g as the two index strings corresponding 
to the two variables contained in the second term of / arc filled in, so that the two variables 
of the first term of g are contained in the second term of /. Adding g to f then has the effect 
of eliminating some disagreements between the two variables as desired, while leaving everything 
unchanged outside the block in question. 

Otherwise, all of the elements marked ? must be 0. We can apply the same argument to the 
terminal string and to the other binomials to obtain g in all cases except where / contains 



1 

1 1 

1 



1 

1 

1 1 



1 1 

1 




1 1 


1 



In this case, because we have a 00 marginal in the first two coordinates of the second term of 
/, we must have one in the first term. If that element contains any 1, by adding a multiple of a 
binomial of degree 2 involving it and the third element in the first term, and then another multiple 
of a binomial of degree 2 involving the third element and the first element in the first term, we can 
construct a g essentially as before which reduces disagreements. The only case where we cannot 
apply this argument to this fourth element is when / contains 



1 


•• 


• 


1 




1 


1 •• 


• 1 


1 


1 


1 •• 


• 1 







1 


•• 


• 








1 •• 


• 1 


1 







•• 


• 


1 





•• 


• 










1 •• 


• 1 






In this case, we let g be this binomial of degree 4, corresponding to switching the middle 
substrings of all O's and all I's. Again, we extend this to g by copying the indices from / outside 
this block to the relevant variables, and we can add this element to / to eliminate this patch of 
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disagreements between the variables in question. This completes the proof by induction. Note that 
we did not use any elements of degree 3 in this process, so Ic„ is in fact generated in degrees 2 and 
4. □ 



This theorem not only shows that the minimal generators are all of degree 2 or 4, but it also 
gives a complete description of these generators. The degree-2 generators come from separations 
of the graph; we will prove a general statement characterizing degree-2 generators of graph ideals 
in Section 5. As for the minimal generators of degree four, we have the following categorization. 

Theorem 3.2. The minimal generators of degree 4 in the graph ideal Ic„ are those elements of 
the form 





1 


B 


1 " 




" Ai 


1 


1 - B 


1 


A2 


1 


1 - B 







A2 


1 


B 





As 





1 - B 


1 




A3 





B 


1 


Ai 





B 


_ 




A, 





1-B 






where the columns correspond to Vi,xi,V2, and X2, Vi and V2 are contiguous blocks of elements, 
and these elements in this order comprise the n-cycle. Here, 1 — B represents the opposite string 
ofB. 

Note that these generators of degree 4 are very similar to the generator of degree 4 in K3, 
namely 



" 1 


1 


1 ' 




" 1 





1 " 


1 










1 


1 











1 







1 


1 


_ 


1 


_ 




_ 





_ 



We will prove a general similarity theorem in this vein in Section 4, when we classify generators 
of a given degree. 

3.2 Complete Bipartite Graphs 

Another nice class of models is the i^m.n model, where Km,n is the complete bipartite graph with 
partite sets of m and n vertices. We first prove the following theorem. 

Theorem 3.3. The graph ideal for G = K2^n is generated in degrees 2 and 4 (for n > 2). 

Proof. As in the proof of Theorem 13. ![ we use binomials of small degree (< 4) to transform a 
binomial in this ideal to one whose two terms share a variable, completing the proof by induction. 
Let the vertices of the two-element partite set heV = {vi, V2}, and let the vertices of the n-element 
partite set he W = {wi, . . . ,Wn}- 

For each monomial M, we define the submonomial Mij, i,j € {0, 1}, to be the product of the 
variables with ij in the index string corresponding to the two-element partite set; we will write 
that index string first. For a monomial M, we define aij{M) to be the total degree of Mij, and 
bij,k,i{M) to be the number of appearances of the digit k in the ^/-position of the index strings 
of the variables in M^j. Here k G {0, 1} and I E {1, . . . ,n}. In other words, these function values 
enumerates the marginal in the direction {vi,V2, wi) with the set values k). We can of course 
recover aij{M) = bij^^i + bij^i^i for any /. 

Then we have the following easy lemma. 
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Lemma 3.4. If Mi and M2 are monomials such that bij^k,i{Mi) = bij^k,i{M2) for all i,j,k,l, then 
their difference can be expressed as a sum of multiples of quadratic elements of Iq- 

We do this simply by, for each Mij , using quadratic generators corresponding to the separation 
{wi, V,W \ {wi}) to move around the bij^i^i{M) I's in the ^th column. Consequently, we need only 
to connect monomials with the same marginals and different -F-values. We introduce an additional 
definition. 

Given a monomial M, the function Cij^i{M) is defined to be the subset of 0, 1 which appears in 
the w;;-position in the variables of Mij. Explicitly, this contains when bij^Q^i{M) > 0, and 1 when 
6ij, i,KM) > 0. 

We now unspool a series of moves designed to connect all of the remaining F-values of monomials 
with the same marginals. 

Lemma 3.5. Suppose we have a monomial M and a column I such that 1 G cqij{M),ciqj{M) and 
£ Coo,i(M), cii^;(M). Then M is equivalent by adding a multiple of a binomial of degree 4 to a 
monomial with the following changes to the aj 's and bj 's: 

+ 1 '■ ^01,0,h ^10,0J) ^00,l,«) 

— 1 : ^oi,i,«) ^io,i,«, ^ii,o,«) ^oo,o,z- 
Proof. This corresponds merely to adding a multiple of the degree 4 binomial 



h 
h 
h 
h 



h 
h 
h 
h 



where the columns are vi,wi, and V2, and all other vertices in some order. This binomial comes 
from the minor given by contracting all the Wi, i 7^ I, into either vi 01 V2. □ 

Lemma 3.6. Suppose we have a monomial M such that for each column I, there exists some index 
il such that ii G coi.i(M), cio,/(M). Then M is equivalent by adding a multiple of a binomial of 
degree 2 to a monomial with the following changes to the aj 's and bj 's: 



+1 : an, aoo> ^oo,i;,/) 



for all I G {1, ... , n}. 



Proof. If / is the index string composed of the i/, this corresponds to adding a multiple of the 
degree 2 binomial 



" 1 


1 


/ ■ 




" 1 


/ " 








/ 







1 / 



where the columns are indexed by vi, V2, and the wi. 



□ 



Now, suppose we have two monomials with the same marginals, that is the same image under 
(j)G- Add multiples of the binomials from Lemma l3.5l and Lemma 13.61 until one can no longer apply 
these; since both increase an + ^ fen,;, one will not go around in circles. Our monomials M and 
are now in "reduced" form, in the sense that neither move can be applied. We break the situation 
down into cases. 
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Case 1. Suppose that Mn and Nn are both not equal to 1, so that both M and N have an entry 
which is 11 in the (^1,^2) direction. 

If, for each there exists an index ii such that ii E cii^i{M) and ii € cii^i{N), then, as desired, 
we can simply extract the variable from both M and N; in other words, for both M and N , 

we can find a monomial with the same values of a and b containing this variable. 

If this is not the case, then there exists an I for which without loss of generality ciij{M) = {1} 
and cii^i{N) = {0}. Looking at the {vi,wi) marginal, there exists at least one marginal 11 because 
of the first condition; consequently, there must exist at least one of these marginals in A^. Since 
cii,i{N) = {0}, the only other option is that 1 € cio,/(A^). Similarly, considering the {v2,'Wi) 
marginal, we must have 1 G cqi^i{N). 

Now, we have ^ii,o,K-^) = k > 0. Since biifi^i{M) = 0, looking at the 10-count in the {vi,wi) 
direction, we must have 6io,o,/(-^) = &io,o,/(-^) + k. However, looking at the 00-count in the {v2, wi) 
direction, it now follows that we must have 600,0,/ (^) = 6oo,o,i(-^^) + ^) ^"^^ particular G cqq^i{N). 
This is a contradiction, since we can now apply a move as in Lemma l!-{.5l to A^, contradicting the 
assumption that N is reduced. 

Case 2. Exactly one of Mu and Nu is equal to 1. 

Suppose without loss of generality that Mu / 1 and A^n = 1. Take any ii G cii,i(M); this ii 
must be in ciq^i{N) and cqi^i{N). This means that we can apply a move as in Lemma 13.61 to 
again contradicting the hypothesis that N is reduced. 

Case 3. Both Mu and Nu are empty. 

In this case, it follows immediately that 6io,i,i(M) = 6io,i,i(iV) for all i and / by considering the 
{vi^wi) marginals equal to (l,i). If Miq ^ 1, this means that we can find an ii for all / such that 
this number is nonzero, and we can then pull the corresponding variable out of both M and N . 
Similarly, if Mqi 7^ 1, we can find a shared variable there, and if both of these are 1 then applying 
the same argument to Mqo = M and A'^qo = N finishes the job. 

Our litany of cases has come to an end, completing the proof of Theorem 13.31 Note again that 
we have not only shown that /i(X2,n) = 4, but also given an explicit generating set in degrees 2 
and 4 for „ • D 

For Km,n where m, n > 2, the answer is less clear. The statement and proof of Theorem 13.31 
indicate that for m fixed, as n gets large, the maximum degree of an element in the Markov basis of 
the graph ideal of i^m,n stabilizes. The degree, on the other hand, certainly goes up as min(m, n) 
does; for instance, there is an element of degree 2m in the Markov basis of Km,m, and we have the 
following result. 

Proposition 3.7. Fix m > 2. Then for n > (^)2™'~^, there is an element of degree 2^~^ in the 
graph ideal of Km,n ■ 

Proof. Let the vertices of Km,n be {vi, . . . , Vm} U {wj}, where I = (ii, . . . , im) is an index string 
of length m, consisting of precisely two I's, and some number of O's and 2's. There are precisely 
g^pj^ strings. 

Specify the marginals as follows: between vj and wj, insist upon ij marginals of 11, 2 — ij 
marginals of 01, 2'""^ — ij marginals of 10, and 2™"^ — 2 + ij marginals of 00. What this means 
is that exactly two variables with coordinate wi equal to 1 occur, and that the sum of the v- 
coordinates (considered as vectors) is precisely I; it furthermore specifies that each of and 1 
occurs 2™~^ times in each Uj-coordinate. 

For each wj, there are only two ways to express the vector / as the sum of two 0-1 vectors. 
Since we consider all index strings /, we obtain that for each diamond in the natural Boolean 
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partial order of binary strings of length m, either the top and bottom elements are in the set of 
u-coordinates of table entries, or the middle two entries are. By an easy induction, it follows that 
the set of ^-coordinates, which numbers only 2"^~^, must consist of either all strings with an even 
number of I's, or all strings with an odd number of I's. From here, we can easily compute the 
t(;-coordinates of each of these entries. 

These resulting tables are the only two which satisfy these marginals, and thus their difference, 
an element of degree 2"^~^, must be in the Markov basis of Ik^^ n ^ desired. For all n > (™)2™~^ 
there is a move of degree 2"^"^ by Corollary 14.21 □ 

We suspect that the following conjecture, an extension of the result for i^2,n7 holds. 

Conjecture 3.8. The graph ideal for G = Km,n is generated in degree at most 2'^™{™-'"} . 



4 Combinatorial Classification of Minimal Generators of Low De- 
gree 

In this section, we give algorithms for computing all generators of a given degree in the graph 
ideal Iq- For degrees two and three, we give an explicit combinatorial characterization of these 
generators, giving a generating set which generates Iq in degree less than or equal to 3; for arbitrary 
degree d, we categorize these generators as pullbacks of a distinguished generator in the graph ideal 
of a fundamental graph X^- The key lemma is the following, relating generators in Iq to generators 
in a minor of G. 

Lemma 4.1. Let G be a graph, and let f = li-Pi^ — Hp/^ be a binomial contained in Iq- Then we 
have the following. 

(a) If Vi corresponds to a column where all the index strings Ij have the same value, then f is 
a minimal generator if and only if f is a minimal generator of the graph ideal G\vi , where f is the 
natural image of f with the column Vi deleted from each index string. 

(b) If Vi and vj are adjacent and correspond to columns where for each index string Ij or Ik, 
the value of that string in each column is identical, then f is a minimal generator if and only if f 
is a minimal generator, where f is the natural image of f with the two columns Vi and Vj fused. 
Here, f is an element of the graph ideal of G with those two vertices fused. 

(c) Suppose Vi and vj are any two vertices and correspond to columns where for each index 
string Ij or Ik, the value of that string in each column is identical. In this case, if f is a minimal 
generator of the graph ideal of G, then f is a minimal generator of the graph ideal of G with those 
two vertices fused. 

Proof. In each case, a decomposition of / into generators of lower degree can be lifted via the 
obvious method to a decomposition of /. In (a), this is simply inserting the shared value of Vi into 
each index string to form a valid index string for G; in (b) and (c), this is simply duplicating the 
value of each index string in the obvious manner. 

In cases (a) and (b), any decomposition of / must necessarily satisfy the property that each 
binomial used has the property in question, by considering in (a) any marginal containing Vi and in 
(b) the marginal corresponding to the edge ViVj. Therefore, a decomposition of / naturally yields 
a decomposition of /, so if / is a minimal generator / must be also. □ 

This lemma has a corollary legitimizing our notion of minor. 
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Corollary 4.2. If H is a minor of G, then fJ,{H) < fJ,{G). 

Proof. It suffices to show IJ-{H) < ^{G) if H is obtained from G by a single vertex deletion or edge 
contraction. However, if it is obtained by a vertex deletion, then by part (a) of Lemma l4. 11 every 
minimal generator of a given degree in Ih lifts to a minimal generator of the same degree in Iq. 
Similarly, part (b) of Lemma l4.ll guarantees that ^{H) < fj,{G) if H is obtained from G via an edge 
contraction. □ 

A natural extension of this is the following, which agrees with the data in Section 2, but which 
we have been unable to prove. 

Conjecture 4.3. If H is obtained from G by deleting an edge, then fJ,{H) < ^{G). 

This set of minimal generators comes with a group action. In particular, the group (Z/2)" 
acts naturally on C[p/], via the element (ci,...,c„) sending a variable Pii-i„ to where 
jj. = + Cr; the sum is evaluated in Z/2. This action consists merely of flipping O's and I's in some 
positions. 

Furthermore, the automorphism group Aut(G) acts naturally on C\pi\ as well, by permuting 
the indices according to the permutation of the vertices of the graph, so we have a natural action 
of (Z/2)"' © Aut(G) on C[pi\. This action maps Iq onto itself; we make the following natural 
definition. 

Definition 4.4. Two generators are equivalent if they lie in the same orbit of C[pi] under the 
action of [TLjlY x Aut(G). 

If two generators of graph ideals and reduce to equivalent generators in a basic graph H 
by means of the above manipulations, we say that they are weakly similar. If they furthermore 
reduce to equivalent generators using only manipulations of type (a) and (b), we say that they 
are strongly similar. We are now prepared to define the object pivotal in our categorization of 
generators of degree d. 

Definition 4.5. Fix a degree d>2. The fundamental graph has vertex set {Si,Ti), where Si 
and Ti are subsets of {1, . . . ,d} with cardinalities \Si\ = |Tj| < d/2, and if \Si\ = d/2 then 1 G Si. 
Two vertices {Si,Ti) and {S2,T2) are connected by an edge if \Si H S2\ = |Ti nT2|. 

To this fundamental graph is associated a distinguished element of Ix^- 

Definition 4.6. Fix a degree d>2. Then the distinguished generator fd £ Ix^ is the binomial 



' h ' 




' Ji ' 


. h _ 







where Iji = 1 if j £ Si and otherwise, and similarly Jji = 1 if j £ Ti and otherwise. 

It is clear that the distinguished generator is actually in Ix^, since for all adjacent Vi and Vj, 
the number of 11-marginals in I is equal to Si R Sj, while the number of 11-marginals in J is equal 
to Ti n Tj. By definition of X^, these are equal, and furthermore, the number of 11-marginals 
determines the numbers of all other marginals (along with the numbers of I's in each column of / 
and J, which are of course identical.) 
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Figure 1: The fundamental graph 



Example 4.7. The fundamental graph X2 has two vertices (1, 1) and (1, 2) which are not connected 
by an edge. The distinguished generator of is the binomial 



" 1 1 ■ 




"10" 







1 



which is just as 2 x 2 determinant. The fundamental graph has nine vertices which are 
(1, 1), (1, 2), . . . , (3, 3). Two vertices (i, j) and (/c, I) are connected if and only if i / /c and j / /. 
Each vertex has degree four: X3 is the edge graph of the 4-polytope A2 x A2, the product to two 
triangles pictured in Figure 1. Note the six triangular prisms which appears as minors of X^. The 
fundamental generator of Ix^ is the binomial 



1 1 1 
1 1 1 
000000111 



1 1 1 
1 1 1 
1 1 1 



With these definitions, we can formulate the main theorem of this section, categorizing all 
generators of degree d. 

Theorem 4.8. Let G he any graph. Then the minimal generators of degree d in Iq can be enu- 
merated by the following procedure. 

1. Consider all graph homomorphisms (j) from minors H of G to X^. 

2. Determine if the natural image of the fundamental generator fd is a minimal generator on 
the image subgraph. 

3. If so, consider the pullback (j)*fd, a binomial of degree d in Iq- We can pull this back to H 
using Lemma \4.1\ part (b) and (c), and then to G using Lemma \4.1\ part (a) and (b). 



15 



A graph homomorphism G ^ H is simply a map (p from vertices of G to vertices of H such 
that if u and v are adjacent in G, (p{u) and (f){v) are adjacent in H. Note that since we have used 
Lemma 14.11 part (c) , the pullback will not always be a minimal generator but the set of moves 
calculated in this way contains all minimal generators Iq of degree d. 

Proof. Suppose we have a generator / of degree d m Iq, written in tableau notation as 



" h ' 




' Jl ' 


. Id _ 




. Jd _ 



where / and J are 0-1 matrices. Each column of I has the same number of I's as the corresponding 
column of J. If this number is either or d, we can delete that vertex via Lemma l4.1l (a) to obtain 
an equivalent generator in a minor of G, proving the theorem by induction. 

Therefore, by flipping with 1 if necessary, we can assume that each column has at most d/2 
I's, and if it has d/2 I's then /j has a 1 in that column. To each vertex Vi associate the pair (S'j, Tj), 
where j S Si if Ij has a 1 in the column corresponding to Vi, and j € Tj if Jj has a 1 in that column. 
If two adjacent vertices have the same pair, we can contract the edge between them via Lemma l4. II 
(b) to again obtain an equivalent generator in a minor of G. 

If this is not the case, then we claim that the map given by sending Vi to the vertex (Si,Ti) in 
Xd is a graph homomorphism. Indeed, all that we need to check is that if Vi and vj are connected 
by an edge, Si n Sj = Tj n Tj (and it is not the case that {Si,Ti) = {Sj,Tj), which is true since no 
two adjacent vertices have the same pair.) But this must be true, since each is just the number of 
11-marginals along the edge ViVj in the corresponding tables / and J, which are equal since / is in 
Ig. 

Furthermore, this map fuses two vertices Vi and Vj only if the corresponding columns of / and of 
J are identical. Therefore, by Lemma l4. II (c). the image of / in //j, where R is the image subgraph 
of Xd, must be a minimal generator of Xj^^. However, this generator is precisely the natural image 
of the distinguished generator in the graph ideal Ir. If this is irreducible, then / is produced by the 
procedure in Theorem 14.81 which we have just done in reverse. If not, then / cannot be a minimal 
generator by the contrapositive of Lemma l4.ll □ 

In this manner, we have reduced the computation of all generators of degree d to the process of 
determining which images of the distinguished generator in subgraphs X^ are minimal generators, 
and of enumerating graph homomorphisms from G to X^i- While the problem of computing graph 
homomorphisms is a difficult one, we can use symmetry techniques to greatly aid us in many cases. 
Theorem 14.81 also divulges which generators are weakly similar: those whose corresponding graph 
homomorphisms have images which are isomorphic subgraphs of X^- 

We now apply Theorem 14.81 to degrees 2 and 3. Consider first generators of degree 2. The 
fundamental graph X2 has two vertices, (1,1) and (1,2), which are not connected by an edge. 
Given a minor of G, it has a homomorphism onto X2 if and only if there are no edges between 
the vertices mapped to (1,1) and the vertices mapped to (1,2). In this case the homomorphisms 
correspond to partitions of the vertices into Vi and V2. Extracting the definition of minor yields 
the following Corollary. 

Corollary 4.9. Let G be a graph. Then equivalence classes of generators of degree 2 of Ig cor- 
respond to partitions Vi U V2 U V3, where there are no edges between Vi and V2. The generator 
corresponding to this is precisely 
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" 101 ■ 




" 111 " 


oil 




001 



where the three columns correspond to Vi,V2, and V3. 

Readers familiar with the study of graphical models and their induced independence statements 
will recognize this as a theorem which says that the only independence statements induced by a 
graphical model are the global independence statements. See, for example fH^ . 

Next, we turn our attention to the case of generators of degree 3, using the methods of The- 
orem ^H] to obtain a combinatorial classification of all such cubic minimal generators of Iq- The 
fundamental graph X3 has nine vertices {(1, 1), . . . , (3, 3)}, with (i, j) connected to {k, I) \i i ^ k 
and j ^ I. By direct computation, the image of the fundamental generator in Ir, R a, subgraph 
of X3, is a minimal generator if and only if R contains one of the six triangular prism subgraphs 
of X3. Therefore, we can classify cubic generators by finding all homomorphisms from minors of 
G to the fundamental graph whose image contains a triangular prism. In particular, if no 
such homomorphism exists, then Iq cannot have a generator of degree 3. This yields the following 
corollary. 

Corollary 4.10. For all n, the graph ideal Ik„ has no minimal generators of degree 3. 

Proof. The only minors of Kn are copies of Km for m < n. Furthermore, all homomorphic images 
of complete graphs are again complete graphs (indeed, of the same degree.) The only complete 
graphs occurring in X3 are K2 and K^, neither of which contains a triangular prism. □ 

Cor ollarv 14 . 1 ( )l shows why edge deletion does not behave well with respect to the graph ideal Iq 
and its Markov basis; this operation obviously can introduce elements of new degrees, since there 
exist graphs with cubic minimal generators, and every graph can be obtained from a complete graph 
by edge deletion. The same technique can be used to show that there are no minimal generators 
of degree 5 in the graph ideal Ik„, since the largest clique in has size 5, and the distinguished 
generator cannot be minimal in the graph ideal of these subgraphs since in fact there are no minimal 
generators of degree 5 in Ki for i < 5. 

In addition to the description of Theorem 14.81 we can obtain a more straightforward combina- 
torial characterization of minimal generators of degree 3 in graph ideals. We start with a pair of 
definitions. 

Definition 4.11. Let G be any graph. Then the 3-coloring graph C3(G) has vertices equal to 
the set of proper 3-colorings ofG. Two 3-colorings are connected by an edge if one can be obtained 
from the other by the following (reversible) procedure: pick a color i € {1,2,3}, pick a connected 
component of the induced subgraph consisting of all vertices with colors not equal to i, and switch 
the other two colors on this component. 

Definition 4.12. A graph G is 3-rigid if the coloring space G-i{G) is disconnected. 

The crucial example of a 3-rigid graph is the same triangular prism that arose in the analysis 
of degree 3 via Theorem 14.81 This graph is 3-rigid, since there are only two proper 3-colorings 
up to permutation of the colors: one triangle {vi.,V2.,v^) is colored (1,2,3), and the other triangle 
(wi, t«2, W'i) must be colored either (2, 3, 1) or (3, 1, 2). Here, the vertex labels are chosen so that Vi 
is adjacent to Wi. It is easy to check that these colorings lie in different connected components of 
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C2,{G)] the connected component of each consists of the colorings obtained from it by permutation 
of the colors. 

We now present a complete description of cubic generators of G based on these combinatorial 
objects. 

Proposition 4.13. Let G he any graph. Then we can enumerate the (equivalence classes of) cubic 
minimal generators of the graph ideal Iq as follows. 

(a) Find all the 3-rigid minors of G. 

(b) For each such 3-rigid minor H , consider all of the connected components of the 3-coloring 
graph of H^. For each connected component, pick a representative 3-coloring Ri. 

(c) For each Ri, i> 1, take the element 





ai 




" &1 " 


f = 


a2 








. "3 . 







where aj^. = 1 if either the vertex k of G is deleted in obtaining H or its image in H is colored with 
color j in Ri , and bjk = 1 if either k is deleted or its image in H is colored with color j in Ri . 

Each of these elements f is a minimal generator of Iq in degree 3, and these elements together 
with the quadratic elements described above generate Iq up to degree 3. 

Proof. Suppose we have a cubic minimal generator of Iq, given by 





ai 




" bi ' 


/ = 


a2 




b2 




. 03 . 




bs 



We construct two 3-colorings Ca and of a minor of G as follows. For the vertex Vj, consider 
the j-th column of A and B. If this consists of all I's or all O's, delete vj. If not, exchange O's and 
I's if necessary so that it has exactly one 1. Then, in Ca (resp. Cb), color Vj with the position in 
which this 1 appears in A (resp. B). 

Next, if two adjacent vertices have the same color in A (equivalent to having the same color in 
B by counting 11-marginals along this edge), then contract that edge. What remains is two proper 
3-colorings of a minor H, and when these 3-colorings are converted to a binomial in Iq via the 
method in part (c), we recover precisely the element /. By Lemma l4.ll the image / in Ih is a 
minimal generator if and only if / is a minimal generator of Iq- 

We claim that the image / is a minimal generator of Ih, i-e. inexpressible as a sum of multiples 
of quadrics, if and only if Ga and Gb are in different connected components of G-i{H). Indeed, 
consider a multiple of a binomial generator, say 





ai 




Cl 


9 = 


a2 




C2 




. "3 . 




. «3 . 



If we convert the two monomials of g into colorings, the same set of vertices will have color 
3 in these colorings. Therefore, the difference between these colorings consists of changing colors 
from 1 to 2 or vice versa. The only way this can be done while preserving the properness of the 
coloring is if the change consists of switching the colors 1 and 2 on some connected components of 
the induced subgraph of H consisting of vertices not colored 3. Therefore, when we add a multiple 
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of a binomial generator, we end 1 up with an element corresponding to a 3-coloring in the same 
component of C^{H), and any two elements connected by an edge in C^{H) differ by a multiple of 
a binomial generator. 

Consequently, the colorings Ca and Cb are in the same component of C^{H) if and only if / is 
not expressible as the sum of multiples of binomials in which is equivalent to / being a minimal 
generator of Ih^ as desired. □ 

Classifying 3-rigid graphs is an interesting problem; the graph C3(G) has been studied in con- 
nection with the problem of picking a random 3-coloring of a graph |16j . Indeed, the flip inter- 
changing two colors on a connected component is precisely the move used in the Wang-Swendsen- 
Kotecky algorithm to pick a random /c-coloring of a graph, and this scheme has ties to mathematical 
physics There are simple operations to produce 3-rigid graphs from other ones, but the trian- 
gular prism seems to be the only 3-rigid graph without a proper 3-rigid minor. 

This method, paralleling Theorem 14.81 can be extended to higher degrees. However, it rapidly 
becomes unwieldy, as the vertices can now be colored with sets of d/2 colors, and the moves are more 
complicated, consisting of all moves keeping one of the colors fixed. For instance, the fundamental 
graph X4 has 34 vertices and understanding the homomorphisms to this graph seems difficult. 



5 Algebraic Degree of Forests 

A recent series of results gave a very thorough description of the family of ideals of decomposable 
models j31ElE]- A special case of these results is the following fundamental theorem. 

Theorem 5.1. The ideal Iq is minimally generated by quadrics if and only if G is a forest. In 
this case, the set of quadratic squarefree binomials in Iq forms a Grobner basis with respect to the 
reverse lexicographic term order with po --o ~< Po-- oi -< ■ ■ ■ ^ Pi-i- 

Sturmfels jTl] posed the natural follow-up problem of calculating the degree of the toric ideal Iq- 
The degree of the toric ideal is interesting in statistics because it provides a natural upper bound 
for the maximum likelihood degree of the toric ideal [E]. In this section we give combinatorial 
formulae for the degree of the graph model ideal Iq when G is a forest. As the maximum likelihood 
degree of a forest is always 1, we see that the degree can be arbitrarily far from the maximum 
likelihood degree. To perform these degree computations, we first recall a result about the degree 
of a general toric ideal. This result can be found in jl3j . 

Theorem 5.2. Let A be a d x n matrix whose toric ideal Ia is homogeneous in the standard Z- 
grading. Then the degree of the ideal Ia (= degree of the projective toric variety defined by I a) is 
equal to the normalized volume of the lattice polytope Q = conv(74). 

Henceforth, the normalized volume of a lattice polytope Q will be denoted Vol{Q). This theorem 
reduces the problem of calculating degree to computing the normalized volume of polytope. We 
now record some some basic facts about the polytope Pq when the underlying graph is a forest. 

Lemma 5.3. The polytope Pq has dimension equal to the sum of the number of vertices and the 
number of edges of the graph (this is true for any graph) . There are precisely 4:- \E{G)\ + 2- \Iso{G)\ 

(j k) 

facets of Pg when G is a forest. If the variables in marginal space are labelled Ui for the variables 

coming from an edge and yf^ for the variables coming from an isolated vertex, then the facets are 
given by the inequalities 
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> and y} ' > 0. 

— Wl; — 



Proof. The dimension formula appears in The polyhedral results are a direct consequence of 
the closed form expressions for maximum likelihood estimates in decomposable models in [TO]. □ 

The following lemma implies that to compute the degree of Iq when G is a forest, one need 
only describe a formula which is valid for trees. Furthermore, this lemma is important for carrying 
out the recursive computations implied by Theorem 15.51 

Lemma 5.4. Let G be a graph with a partition of the vertices {Vi,V2} such that there is no edge 
in G between the Vi and V2- Let Gi and G2 be the corresponding induced subgraphs. Let di and d2 
be the corresponding dimensions of the polytopes Pqi and -PG2/ ^-S; di = \Vi\ + \E{Gi)\ . Then 
we have 



Proof. With these restrictions on the graph G, we have Pq = Pgi ^ Pg2- Equation 15.41 is then the 
usual formula for the normalized volume of the direct product in terms of the normalized volumes 



We now come to the main theorem of this section. 

Theorem 5.5. Let G be a tree with n-vertices. Then the degree of the toric ideal Iq can be 
calculated by the formula 



where the notation G\e denotes the graph G with the edge e removed. 

Proof. As previously indicated, we prove the theorem by calculating the volume of the correspond- 
ing polytope Pq. Theorem 15.11 implies that the pulling triangulation of Pq induced by the reverse 
lexicographic term order above is unimodular. This in turn, implies that the normalized volume 
of Pg is equal to the sum of the normalized volumes of the facets of Pq which do not contain the 
vertex indexed by the variable pi-.-i (see ^Hl for all the definitions and relevant theory). This is 
where the polyhedral description of Pq when G is a forest becomes essential. We see from the 
polyhedral characterization that there are exactly n — 1 facets of Pq which do not contain this 
"last" vertex, and that they are indexed by the edges of G. We will show that the normalized 
volume of the facet of Pg which is indexed by the edge e has volume precisely ^VoI{Pq\^^) which 
will complete the proof. 

There are two cases to consider: either the edge in question has one node a leaf or not (the case 
of the graph on two vertices with a lone edge is clear, by a direct calculation). We will handle the 
two cases separately. 

Case 1: We may suppose without loss of generality that our edge e is {1,2}, the vertex 1 is a 
leaf and the vertex 2 has the edge {2,3} incident to it. Then the matrix whose columns correspond 
to the vertices of Pq has the top eight rows which look like 




of the components of the product. 



□ 
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with this 8x8 block repeated 2" ^ times across the first eight rows. We claim that the facet which 

(1 2) 

corresponds to the inequahty yll > (corresponding to the fourth row of the above matrix) has 
volume equal to ^VoI{Pg\{i^2})- Note the the vertices of Pq which lie on this facet are precisely the 
3 • 2"~^ vertices of Pq which have a zero in the fourth row. First we show that this facet is naturally 
isomorphic to a sub-configuration of -Pg\{i,2}- Consider the matrix ^g\{i,2} whose columns give 
the vertices of Pg\{i,2}- This matrix has two fewer rows than the matrix Aq above and is almost 
the same: its first six rows look like 
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To show the natural isomorphism mentioned above, it suffices to show that there is a unimodular 
linear transformation from the first six columns of the first matrix above to the first six columns 
of the second matrix above. Such a linear transformation is obtained by replacing the first row by 
the sum of the first and second rows, and then deleting the second and fourth rows. We can delete 
the second and fourth rows because they are linear combinations of other rows and hence do not 
change the polyhedral description. Such a transformation is clearly unimodular. 

Now that we have shown that our configuration of 3 • 2"~^ points sits naturally inside -Pg\{i,2}! 
we wish to compute its volume. For this, we show that there is a hyperplane which divides -Pg\{i,2} 



into two congruent pieces, one of which is the convex hull of our new configuration of 3 • 2*^ 
This hyperplane is given by the equation 



pomts. 



(1) (1) , (2,3) , (2,3) (2,3) (2,3) 

Vo -vi' + vofl + 2/0,1 - yi,o - yi,i 



0. 



Note that exactly 2"~^ vertices of -Pg\{i,2} li^ on this hyperplane (these are the ones corresponding 
to the middle four columns of the submatrix of Ag\{i,2} displayed above) and the remaining 2"^^ 
vertices are split equally on each side of the hyperplane. In particular, all of the vertices from our 
configuration of 3 • 2""^ points lie on the nonnegative side of this hyperplane and the remaining 
2"~^ points are on the negative side. Furthermore, there is a natural reflexive symmetry across 
this hyperplane. To complete the proof, we must show more: not only is the point configuration 
naturally "cut in half" by this hyperplane, but so is the polytope -Pg\{i,2}- This follows from a 
direct computation with the eight points listed above. We performed the computation using the 
program PORTA [T]. 

Case 2: The argument is very similar to case 1; we will sketch the relevant details. We may 
assume that our edge is the edge {2,3}. Since neither 2 nor 3 is a leaf we may assume G also 
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contains the edges {1,2} and {3,4}. With these conditions, the first 12 rows of our matrix looks 
like 
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with this block repeated 2"~^ times. We wish to show that the facet defined by the inequality 
y(2^3) y g (corresponding to the eighth row in the above matrix) has volume equal to ^VoI{Pg\{2,3})- 
The vertices which lie on this facet are precisely the 3 • 2^~'^ vertices with a zero in the eighth row in 
this matrix representation. First we show that this facet naturally appears as a sub-configuration of 
Pg\{2,3}- This can be achieved by applying a unimodular transformation to the configuration: the 
key point is that once we restrict attention to the facet, the middle 4 rows of the above configuration 
can be written as linear combinations of the other rows and hence are redundant in terms of the 
polyhedral description. 

Now we show that there is a hyperplane which divides the polytope -Pg\(2,3) iii half. This is just 
the hyperplane given by 



(1-2) (1,2) , 
yo,0 - VO,! + ^1,0 



(1.2) (1,2) , (3,4) , vu,. 
- - yi,i + 2/0,0 + ^0,1 



,(3,4) _ (3,4) _ (3,4) _ 

0,1 yi,o ^1,1 ~ 



Note that our configuration of 3 • 2""^ points are precisely the points on the nonnegative side of this 
hyperplane. Furthermore, there is a natural reflexive symmetry across the hyperplane. A direct 
calculation shows that this hyperplane not only separates the point configuration, but also divides 
the polytope into two symmetric pieces with the the desired integral points as vertices. Thus we 
deduce the desired equation of volumes. □ 

For some special families of trees we use this recurrence relation to deduce simple formulae for 
the degree. These appear in the following corollaries. 



Corollary 5.6. Let Ki^n denote a star graph with n leaves. Then deg{lKi , 



(n!) 



n2_ 



Proof. Removing any edge of the graph Ki^n produces the graph consisting of the disjoint union of 
a Ki^n-i and an isolated point. Hence from theorem 15 . 51 and lemma 1231 we deduce that 



1 /2n^ 

deg{lK,„) =n - - [ ^ ) • deg{lKi, 



r? ■ deg{lK,n-i)- 



Since deg{lKi = 1 we have the desired result. 



□ 
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Corollary 5.7. Let T„ denote the graph of the n-chain and dn = deg{lT„)- Then dn satisfies the 
recurrence 




with di = l. Furthermore, we have the equality of generating functions 

oo , 

^(2n-l)! ^V2^ 

Proof. The recurrence relation 15.71 follows immediately from the formula in Theorem 15.51 and by 
applying Lemma 15.41 to the graph consisting of two disjoint chains of length i and n+1 — i. 

To deduce the equality of generating functions, let y = ^ Yl^=i (2n-i)! ^^"~^" recurrence 
relation implies that 2y' — 2 = y^. Solving the differential equation yields the desired formula. □ 

The recurrence relation and generating function in the case of the n-chain also appears in a 
paper by Poupard JJ, but we do not know how to show that the objects we are counting (simplices 
in a regular unimodular triangulation) are in bijection with the objects she was counting (types of 
binary trees). 
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